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Abstract 

We review different definitions of the current density for quantized 
fermions in the presence of an external electromagnetic field. Sev¬ 
eral deficiencies in the popular prescription due to Schwinger and the 
mode sum formula for static external potentials are pointed out. We 
argue that Dirac’s method, which is the analog of the Hadamard point¬ 
splitting employed in quantum field theory in curved space-times, is 
conceptually the most satisfactory. As a concrete example, we discuss 
vacuum polarization and the stress-energy tensor for massless fermions 
in 1-f 1 dimension. Also a general formula for the vacuum polarization 
in static external potentials in 3-1-1 dimensions is derived. 

Keywords: Quantum electrodynamics in external potentials; current 
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1 Introduction 

Quantum field theory in external potentials has received increasing interest 
in recent years, triggered in particular by the envisioned high intensity laser 
sources, cf. [T] for a review, and the new facilities for heavy ion collisions, in 
which one might probe the chiral magnetic effect [2]. Further applications 
are vacuum polarization corrections to atomic energy levels, cf. the review 
[3] . The topic has a long history, dating back to the seminal works of Dirac 
[3| and Heisenberg and Euler [5]. Schwinger basically brought the theory to 
its present form, by introducing the proper time formalism and the concept 
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of the effective action [B]. For textbooks and reviews on the subject, we 
refer to [ZlIHlEllinilll]. 

There are two types of observables in quantum electrodynamics in ex¬ 
ternal potential that one might be interested in, namely global and local 
ones. The prime example for a global observable is the expectation number 
of electron positron pairs produced by the external field. An example for a 
local observable would be vacuum polarization, i.e., the expectation value 
of the charge density in a static external potential. More generally, the local 
current density j^{x) should be an observable. 

In principle, both types of observables can be studied in the framework 
of effective actions. The particle production rate is related to the imaginary 
part of the effective action, whereas the expectation value of the current 
can be obtained by variation w.r.t. the external potential A^. However, 
there is the practical difficulty that the effective action can be explictly 
computed only for very special forms of the external potential: To the best 
of our knowledge, the most general prescription to compute the effective 
action was given in |12] . where it is implicitly assumed that the external 
potential is analytic. But an analytic potential can not be varied locally, so 
the expectation value of the local current density can not be computed in 
that way. 

On the conceptual side, it would be desirable to have a definition of the 
current density that is independent of the choice of a state. This would 
allow to evaluate it also in thermal states, or in situations where there is no 
preferred vacuum state, due to non-stationarity of the external fields in the 
asymptotic past (e.g., as for a plane wave). 

Classically, the current density of a Dirac field is given by 

j^{x) = —e'4){x)^^'ip{x). ( 1 ) 

Upon quantization, one meets the immediate problem that point-wise prod¬ 
ucts of quantum fields are not well-defined. In the absence of external poten¬ 
tials, one typically proceeds via point-splitting w.r.t. the vacuum two-point 
function, i.e., normal ordering. However, it is not clear how to replace the 
vacuum two-point function in the case of non-trivial external potentials. 

The first proposal for a point-splitting definition of the current seems to 
be due to Dirac [1] . He assumed that in a physically reasonable state w the 
two-point function 

uj{ip{x)ijj{y)) 

has the same singularities as the vacuum two-point function, but possibly 
with modified coefficients. The proposal for a renormalization prescription 
was then to subtract the singular part. As discussed below, this is, in mod¬ 
ern terminology, the assumption that the two-point function is of Hadamard 
form and that the correct renormalization prescription is by point-splitting 
w.r.t. the Hadamard parametrix. This is the renormalization prescription 
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that we are going to employ here. The main shortcomings of Dirac’s treat¬ 
ment are the missing discussion of the inherent ambiguities and the only 
implicit determination of the coefficients of the singular part of the two- 
point function. Both were treated subsequently by Heisenberg [13], who 
showed that there is an ambiguity in the definition of the current, namely, 
one may add multiples of the external current (which corresponds to a charge 
renormalization). Dirac’s method was also used by Heisenberg and Euler [5|, 
Uehling [T3|, and Berber m in their seminal works on vacuum polarization. 
However, it seems to have been essentially abandoned afterwards, being su¬ 
perseded by Schwinger’s prescription. The reason for this may have been 
practical: The singular part is explicitly given in position space, whereas 
the two-point function is typically expressed as an integral over modes, from 
which the smooth remainder is difficult to extract, at least beyond the ap¬ 
proximation of first order in the external potential^ One of the results of 
the present paper is a method, applicable in the case of static external po¬ 
tentials, to subtract the singular part directly in the mode integral, thus 
allowing for more reliable numerical calculations. 

Another dehnition of the local current was given by Schwinger [6|, who 
proposed to perform the limit of coinciding points in Q in a symmetric way, 
leading to cancellations of divergent parts. However, not all divergences are 
canceled in this way: As we will see, a logarithmic divergence remains, un¬ 
less the external current vanishes. Furthermore, even if the external current 
vanishes, Schwinger’s prescription in general depends on the way the limit 
of coinciding points is performed. It should also be noted that some au¬ 
thors seem to omit the parallel transport that is implicit in Schwinger’s 
point-splitting procedure, cf. the footnote in m- Then, further divergences 
remain, and the result may fail to be gauge invariant. It should also be 
noted that, contrary to Dirac’s method, Schwinger’s prescription can not be 
used to renormalize the stress-energy tensor. 

Finally, for the case of a static external potential, there is a dehnition 
for the vacuum polarization based on the mode expansion: 

. ( 2 ) 

This expression seems to originate in the work of Wichmann and Kroll m 
and is widely used in textbooks and reviews, e.g., mm- Here the hrst sum 
is over positive and the second sum over the negative frequency modes, with 
the V’n’s properly normalized. Unless the system is spatially conhned, the 
sums have of course to be replaced by integrals. Furthermore, the individual 

^This is also evidenced by the fact that the actual calculations of Heisenberg and Euler 
|5] proceed by a cut-off in momentum space, not by a point-splitting. They claim to also 
have done the calculation via point-splitting, but remark that it “turned out to be very 
complicated.” 
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sums are divergent. In typical situations, positive and negative frequency 
modes can be paired such that the remaining sum is convergent (this was 
done in [TH] for the case of a charged scalar field). As we will show by an 
elementary example, this expression in general leads to incorrect results. 
The same is true for an analogous expression for the total energy. 

In the context of quantum field theory on curved spacetimes, one also 
encounters the problem of defining Wick powers, i.e., non-linear local oberv- 
ables, in the absence of a preferred vacuum state. In recent years, there has 
been tremendous progress in that field, triggered by the application of the 
tools of microlocal analysis [20] and by the introduction of the framework 
of locally covariant field theory [111122]. In the latter, one does not define a 
quantum field theory on one particular spacetime, but on all possible ones 
simultaneously, in a coherent way. The same applies to observables, such 
as the stress-energy tensor. This is non-trivial, as it necessitates a point¬ 
splitting, which is required to depend only on the local geometric data, in or¬ 
der to fulfill the condition of local covariance. The framework was applied to 
conceptual as well as practical problems, such as the proof of spin-statistics 
|23j and CPT |24j theorems, the discussion of quantum energy inequalities 
[25], or to cosmological back-reaction effects [26]. For a recent review, we 
refer to [27] . 

The framework of locally covariant field theory can be straightforwardly 
generalized to fields that are charged under a gauge group G, in the pres¬ 
ence of an external potential [2S] . Specializing to Minkowski spacetime and 
G = U{1), one obtains a framework for quantum field theory in external 
electromagnetic potentials, in which all observables are constructed in a 
locally gauge covariant manner. Considering the electromagnetic field as 
purely external, one is dealing with a free quantum theory, in which only 
Wick powers need to be renormalized. This proceeds via point-splitting 
w.r.t. the Hadamard parametrix, i.e., as in Dirac’s prescription. We note 
that the well-known chiral anomalies can be understood (and also computed 
in a rather elementary way) as due to the Hadamard parametrix being a bi¬ 
solution to the equation of motion only up to smooth remainders |29j . This 
may be seen as another indication of the appropriateness of the framework. 

One aim of this paper is to review Dirac’s point-splitting procedure in 
modern terms and to relate it to Schwinger’s prescription and the mode 
sum formula (§. We also discuss the remaining ambiguities. This is done 
in the next section. In Section [^ we apply the procedure in a case that is 
analytically tractable, a massless fermion in two space-time dimensions in a 
constant electric field and confined to a spatial interval. This example also 
nicely illustrates the failure of the mode sum formula (§. For applications in 
four space-time dimensions, one meets the inconvenience that the Hadamard 
parametrix is known in position space, whereas the two-point functions of 
interesting states are typically given by mode integrals, where it is difficult to 
subtract the position space singularity. In Section [^ we propose a solution 
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to this problem for the case of static external potentials, and exemplify it by 
considering the Coulomb potential. This method can also be modified for 
the case of a homogeneous time-dependent electric field, which is discussed 
in a separate publication m- 

Notations and conventions 

We work with signature The electric charge is denoted by e 

and the covariant derivative is given by + ieA^. Throughout, we 

work in rationalized natural units, i.e., c = h = 1 and With 

the Lagrangian 

L = - m)^lj - 

these conventions determine the sign in Q. When considering functions 
(or sections) in two space-time variables, the second one is often denoted 
by x'. Hence, a primed covariant derivative is implied to act on the second 
variable. The coinciding point limit of such a smooth function is denoted 
by square brackets, [f]{x) = f{x,x). The spatial part of a Lorentz vector x 
is denoted by x. 


2 Hadamard two-point functions and the parametrix 

In the following, we assume, as Dirac |5, that the physically relevant states 
have two-point functions of Hadamard form, i.e., they are, in four space-time 
dimensions, locally given by 




Wi^{x,x') 


+ 


W2^{x,x') 


+ W 3 ^{x,x)\og 


{{x — (x — x')| 


— (x — x') 


A2 




+ w^i{x,x'), (3) 


with A and B spinor and co-spinor indices, 

(x — x')1 = (x — x')^ — ie(x — x')°, (4) 

A a length scale, and smooth functions Wi. The sum of the first three 
terms, i.e., H, is usually called the Hadamard parametrix. Before discussing 
it further, let us recall the motivation for the assumption that the two-point 
function is of the above form. 

Dirac’s motivation was that the vacuum two-point function is of this 
form, and that the introduction of an external potential should not change 
the structure of the divergences. Nowadays, one can give further motivation, 
in particular due to the lessons learned in QFT on curved space-times. To 
begin with, the Hadamard form propagates, i.e., if the two-point function is 
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of Hadamard form in the neighborhood of some time-slice and the potential 
is smooth, then it is of Hadamard form everywhere^ Hence, for a potential 
that is switched on at some finite time, the two-point function that coincides 
with the vacuum two-point function in the past is of Hadamard form. In par¬ 
ticular, this can be used to show that a two-point function of Hadamard form 
always exists [28], by adapting the deformation argument of |34) . Another 
relevant result is that, on Minkowski space and for generic smooth static 
external potentials, the ground state is of Hadamard form |35|. Finally, 
in states with two-point function of Hadamard form and smooth truncated 
n-point functions, the higher moments of smeared Wick powers are finite 
|36j In particular this means that the Hadamard condition ensures finite 
fluctuations of the smeared current density. 

Let us now review the Hadamard form of the two-point function. As 
suggested by Radzikowski [20] , it can be characterized by the condition^ 

{x)iPa{x')) + = iS\{x, x'), (5) 

( 6 ) 

and the requirement that the two-point function is of positive (negative) 
frequency in its first (second) argument. The latter is made precise by a 
condition on the wave front set (for an introduction to this concept we refer 
to [OH])- In Q, S' = S®''^'" — S^'®* is the difference of advanced and retarded 
propagator of the Dirac operator (including the external potential), and in 
(§, u and V denote test sections of the spinor and the co-spinor bundle, 
with which ip and ijj are integrated. 

One can show |33l [28] , that these requirements fix the two-point function 
up to a smooth part, denoted by W 4 in Q. Furthermore, the first three 
terms in Q, are constructed locally and covariantly out of the geometric 
data, or, more precisely, the coinciding point limits 

[VMi---V;.,v;^...v;iy,], (7) 

of covariant derivatives of the Wj, 1 < f < 3, are tensors constructed out of 
the external field strength F. 

The proposal of Dirac [1] , which was later “rediscovered” in the context 
of QFT on curved spacetimes [39llini[lllll2[|2T] J^is then to define non-linear 
local observables (Wick powers) by a point-splitting in which the Hadamard 

^The first such result was obtained in m for scalar fields and generalized to Dirac 
fields in [3^. A modern proof using the tools of microlocal analysis was given in [^. It 
is straightforwardly generalizable to the case of non-trivial external potentials, 
partial converse of this statement was recently proved I37|. 

"^A further condition is that the two-point function is the two-point function of a state, 
i.e., positive, ui{ip{u)'il>{u)) > 0 for u a test section of the spinor bundle. 

®Marecki, apparently unaware of Dirac’s work, but inspired by QFT on curved space- 
time, “reintroduced” this into QED in external potentials |43 |. 
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parametrix H, i.e., the first three terms in ([^, are subtracted. Concretely, 
for the evaluation of the Wick power in the state uj, one obtains 

In particular, this provides a definition of the current j^, by tracing the 
above with 7 ^. It is clear that for the computation of the coinciding point 
limit of W 4 ^ only a finite number of the tensors 0 have to be known (up to 
four derivatives for Wi, up to two derivatives for W 2 and no derivative for 
Vks). If the Wick power contains derivatives (as needed for the stress-energy 
tensor), more of these tensors have to be known, but always just a hnite 
number for any fixed order of derivatives. 

Before we finally compute the tensors Q that are relevant for the def¬ 
inition of the current, let us briefly comment on the ambiguities in the 
definition of Wick powers. These were not treated by Dirac, but a detailed 
discussion in the context of scalar fields on curved space-time can be found 
in [ 211111 ]. One obvious ambiguity is the choice of the length scale A in 
the logarithmic term in the parametrix (a change would amount to adding 
a multiple of \W‘i^{x) to the Wick power Tp a{x)) ■ More generally, one 
may allow for the freedom to add to a Wick square any local tensor of the 
correct dimensionj^ For the Wick square in four space-time dimensions, 
these are multiples of m? and R. For the Wick square ip^pJA on Minkowski 
space, examples of such tensors would be or This 

freedom, however, is constrained by the requirement that the current should 
be conserved. Once the Wick square is chosen such that the current is 
conserved, the only remaining ambiguity consists of adding multiples of the 
external current [28]. This corresponds to a charge renormal¬ 

ization. Furthermore, as we will see below, in four space-time dimensions, 
this can be implemented by a change of the length scale A. 

2.1 The explicit form of the parametrix 

We want to explicitly compute the tensors 0 that are relevant for the point¬ 
splitting renormalization of the current density. Readers not interested in 
the details of this calculation, may skip towards the end of the subsection. 

Due to the condition Q, it is advantageous to consider the singular 

behavior of the retarded and advanced propagators of the Dirac 

operator — m. To obtain these, we follow the strategy of Dimock |45| : 
The operator 

P = — m) = -|- -|- 

is normally hyperbolic, i.e., its principal symbol is the metric and there are 

no first order derivatives. For such operators, there are unique retarded and 

®This can be seen as adding to the Hadamard parametrix H a smooth term that is 
locally and covariantly constructed out of the geometric data. 
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advanced propagators, which are formally given by [46j 


^ret/advj-^^ ®') = ^ 14 (x, {x, x'), (8) 

k=0 

where the Hadamard coefficients 14 are smooth and the Riesz distributions 
Rj in four space-time dimensions are given by 

^ret/adv(^^^,) ^ ^0(^±(^x - x'f)6{{x - x'f) , 

Rl%t^ix,x') = 22 /c+i^fc!(fc_i), (a^ - x'f^^-^'^9{±{x - x'f)e{{x - x')^), 

where in the second line k > 1. The Hadamard coefficients can be deter¬ 
mined recursively by solving the transport equation 

(x - x'yV^Vk{x,x') + kVk{x,x') = -kPVk-i{x,x'), (9) 

with the initial condition Vo(x,x) = 1. This proceeds by integration of 
PVk-i along the straight path from x to x^ For example, one straightfor¬ 
wardly obtains that 

Vq{x,x') = eyi'p (^ie J Afj_{x'+ t{x — x')){x — x')^dt 

is the parallel transport along the straight path from x' to x. For the deter¬ 
mination of the tensors Wi, 1 < i < 3, we will need the coinciding point limit 
of up to three covariant derivatives of Vq and up to one covariant derivative 


of 14 . For the parallel transport, one easily obtains 

N = 1 , ( 10 ) 

[V;.Ho] = 0, (11) 

[V;,V,Ho] = (12) 

[V^V^VaHo] = fe {d^F,x + d,F^x). (13) 

For Vi, we obtain, using (|^ and (10) - (13), 

[l"i] = -[PVo] = -ie^YF^, - m\ (14) 


= -W^uPyo] = -le^^Yd^Fxp - ied^Fpx. (15) 

It is clear that in principle one can compute coinciding point limits of arbi¬ 
trarily many covariant derivatives of I 4 ia this way. 

The expansion ([^ is formal in the sense that the series does in general not 
converge; however, it is asymptotic in the sense that if it is truncated after 
k = N, then the difference of the actual and the truncated 

propagator is of regularity and can be bounded by |(x — x')^i'^“^i| 
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j46| . As we are ultimately interested in the limit of coinciding points, the 
formal expansion Q is completely sufficient for our purposes, also for the 
similarly defined parametrix. 

The retarded/advanced propagator for the Dirac operator — m\s now 
defined as 

^ret/adv ^ 

It is obvious that these have the required support property and that 

= 6 . 


That also 

^ret/adv o - m) = 6 

was shown by Dimock [35]. In particular, it follows that 

(17) 

where the operator is the conjugate Dirac operator (for charge —e) 
acting from the right on the second variable (the change of sign is due to 
partial integration), i.e., 

f*'F{x, x') = (d'^ - ieA^{x'))F{x, x')-i^. 


Knowing the form of the retarded/advanced propagator for the Dirac 
operator, we may now tackle the singular structure of the two-point func¬ 
tion, i.e., the Hadamard parametrix. The crucial point is that there are 
distributions fulfilling 

T+ - T- = 2TTi{Rf'' - Rf) (18) 

for j = 2 + 2k, k G No, and which are of positive/negative frequency in their 
first argument (again, this can be made precise using the concept of the 
wave front set). Explicitly, in four space-time dimensions, they are given by 






-1 


2vr (x-x')L’ 


-1 


22fc+i7rA:!(fe-l)!^ ’ ^ PF 


where A is a length scale and we used the notation Q. Note that 
follows straightforwardly from 


. = P-F i-n5{x), 

X±l£ X 


log(x ± ie) = log |x| ± iTT9{—x). 


(19) 

( 20 ) 
PI 


Hence, we may now define the parametrix h for the operator P as 


h^{x,x') = —^Vkix,x')T^_^^f^{x,x'). 
^ k=0 
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Analogously to the construction (16) for the retarded/advanced propaga¬ 
tor for the Dirac operator, we may define the parametrix H for the Dirac 
operator as 


H^{x, x') = I (^-if + if*' - 2m) h^{x, x') (21) 


Here we took the mean of the action of the auxiliary Dirac operator —i 
on the left and the right in order to fulfill ([^J^ As a straightforward 
quence of (18), we then have 


f — m 
conse- 


H^{x,x')—H {x^x) = iS{x^x). 


( 22 ) 


By construction, is now a solution modulo smooth terms to the 
Dirac operator if — m in the first variable, and to its transpose in the 
second variable, fulfilling ([^ and the positive frequency condition. By the 
uniqueness result mentioned before 0- it follows that for any state cu with 
Hadamard two-point function, we have 

uj{Tp^ {x)ipA{y)) = H+{x,y)^ + R^{x,y), (23) 

uj{ijjA{y)'ip^{x)) = -H-{x, y)^ - R^{x, y), (24) 


where R is smooth. That the two smooth terms on the r.h.s. of these equa¬ 
tions add up to zero is a consequence of ([^ and (22). 

We are now finally in a position to state the functions Wj, 1 < i < 3 
explicitly: 


Wi{x,x') = -^7^(x - x'YVq{x,x'), 

W 2 {x,x') = ^{if - if*' + 2m) Vo{x,x') + ^{x - x'Y{-i^,Vi{x,x')], 
Wz{x,x') = ^ (^f -if*' + 2m^ Vi{x,x') + 0{x-x'). 

From these, we obtain the relevant coinciding point limits of their deriva- 

^We recall that the action of the auxiliary Dirac operator on the retarded/advanced 
propagator gj^gg ^ret/adv^ regardless of whether one acts from left or right, cf. 

0 - 
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tives: 


[Wi] 

= 0, 

1 



[V^lTi] 

= - 

i 



[V(^V,)Wi] 

= 0; 

1 



[V(^V,V;,)Wl] 

= 0, 

1 



[V(^V,VaV^)Wi] 

= 0, 

1 



[ 1 C 2 ] 

= - 

1 

477^ 

m, 



_ i 

1677^ 


■{'yv^'y^Y}Fxp - ‘i-lprr?) 

[V(^V,)W2] 

= - 

e 

2477 




+ 

e 

3^ 



[VDs] 

= - 

1 

1077 

- 

- ^em-t^YFxp - 4g^ 


Note that for the last equation we used Synge’s rule, i.e., 

[v;,Fi] = V,[V^] - [V,Vi]. 

These coincide with the results of Heisenberg [13] (note, however, that there 
the terms with more than one 7 matrix are omitted, as they do not contribute 
to the current). When one is interested in the current these have to be 


traced with 7 ^. The non-zero traces are 

tr 7 ^[V^iTi] = (25) 

tr 7 ^[V^iT 2 ] = (26) 

tr 7 ^[V(^V.)W 2 ] = (27) 

tr7^[W3] = -i^9"T^A. (28) 


2.2 Comparison with Schwinger’s prescription 

With our explicit knowledge of the form of the two-point function, we may 
now compare Dirac’s method to Schwinger’s | 6 | prescription. The latter 
consists in computing 


j^(x) = I lim (yo{x,x + + tv)ipAix) 

-Vo{x -b tv, x)'4)a{x + tv)i)^{x)) 7 ^ 5 , (29) 

where n is a future-pointing time-like unit vector. Before considering this 
any further, let us pause to remark that the parallel transport Vq was not 
explicitly included in Schwinger’s original work | 6 |. Only the later work 
|16j states that it has to be contained. Nevertheless, the parallel transport 
seems usually not to be taken into account, cf. for example the textbook and 
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review mM- This leads to a breakdown of gauge invariance and spurious 
divergences!^ 


Let us now evaluate (29) in a Hadamard state a;. Using (25) - (28), we 


see that the most divergent expressions cancel, and we are left with 


lirnlogt 


t^O 
2 


+ v'^dxFp^ - ^vpv dPF^p, 


where R was the smooth remainder in (23), (24). The first term on the 


r.h.s. is the one that Dirac’s method gives. The second term is a logarithmic 
divergence, which is not canceled, unless the external current vanishes. Then 
also the last term vanishes. Even then, the third term remains. It shows 
that the definition depends on the direction in which the limit is performed^ 
Its presence also spoils the conservation of the current. To eliminate this 
term, one would have to consider a different limit for each of the components 
of jp, by choosing v to be the unit vector in /r direction for the computation 
of the ^ component (the cancellation does not depend on the fact that v 
was originally constrained to be time-like). 

The cancellation of the most severe divergences rests on the fact that 
these are anti-symmetric. However, this is only the case if the parallel 


transport is taken into account, as in (29). Otherwise, spurious divergences 
remain, e.g., for the leading divergence, 


vP'jpVoix + tv,x) 

F 


vP'yp,Vo{x,x + tv) 

F 


(iv'^eAy 


i[v 


'eA„ 


+ 0{t) 


(30) 


where for simplicity we assumed that A^ is constant in u-direction. There 
is also a finite remainder, which is obviously not gauge invariant. The occur¬ 
rence of such spurious divergences seems to be the reason why for example 
in m more than one renormalization condition for the current is imposed 
(even though we know that there is only one ambiguity). 

®Let us also mention that the textbook [7] contains, in Section 9.7, a further mistake, 
stating that the above symmetric time-ordered limit should be performed as 

I lirn ^lirn^ [V’(®° + t,x + s)'tp{x°,x) — ijj{xP,x)ijj{xP — t, a; -|- s)] , 

where s is a spatial separation. But after taking the limit f —>■ 0 for Hxed s 7^ 0, the fields 
anti-commute, yielding 

lim ip(x°,x + 8)ip(x°,x), 

which is divergent. 

®For a coinciding point limit from a lightlike direction, this term would actually diverge. 
Analogous direction dependent terms were also encountered in ES] in the renormalization 
of the stress-energy tensor. 
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Finally, let us investigate the mode sum formula ([^ for the vacuum 
polarization. Assume a static external electric field, in static gauge Aj = 0. 
Let be the normalized modes of the Hamiltonian 

H = —+ m + 7°eAo 

and ipk = Here kis a generalized momentum labeling the eigenstates 

(spin degrees of freedom are also included). Then the vacuum two-point 
function is given by 





V'f (^)V’A,fc(y)e (31) 

V’f (32) 


where x = (x°,x). Choosing v = e^, the unit vector in time direction and 
computing the vacuum polarization according to (29), we obtain 


^(Pix)) = 7 lim 

2 t->0+ 



If it were allowed to interchange the integral and the limit here, one would 
obtain the mode sum formula ([^. However, the integrals are not absolutely 
convergent, so this is not possible. In some situations (we will encounter 
one below), the positive/negative energy modes for a given generalized mo¬ 
mentum k have exactly opposite energy and the same modulus IV’fcl- The 
mode sum formula (§ would then yield a vanishing vacuum polarization. 
The above expression, however, gives 


ie lim sin(teAo(x)) / \ipk{x)fe^^^^d^k, 

JEk<0 ^ 

which does in general not vanish, as the integral diverges as t —>■ 0. We will 
discuss an explicit two-dimensional example below. This again shows that 
neglecting the parallel transport in Schwinger’s prescription in general leads 
to wrong results. 


3 An explicit example in 1+1 dimensions 

We consider the massless Dirac field in I -|- I dimensions, confined to the 
spatial interval [—L/2,L/2] in the presence of a constant electric field E. 
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We choose static gauge, Aq = Ex^, Ai = 0, and the 7 matrices 


7° = 


0 r 

1 0/ ’ 


7 ^ = 


0 1 
-1 0 


The Dirac equation 


may be written as 


if'll; = + ieAf)'il^ = 0 


ido4) = i 


. f-l 0 


0 1 


diij) + eEx^4>, 


(33) 


with (j) = 

In order to construct the two-point function, we proceed by the mode 
decomposition. The r.h.s. of (33) is written as Hp in the following. Eigen¬ 
vectors of H for eigenvalue k are 


= ( 0 ) ^ 




It follows that solutions to if ip = 0 are given by 

^-^eE{x^P+ik{x^+x°) ^ ^ ^^eE{x^)^+ik{-x^+x°) ^ 

To ensure self-adjointness, we have to impose boundary conditions. It does 
not seem to be physically reasonable to impose (anti-) periodic boundary 
conditions, as a particle would be accelerated around the loop, taking up 
energy each time. It seems more appropriate to impose bag boundary con¬ 
ditions m that ensure the vanishing of the current at the boundary 

nVliL/2 = ±^/’l±L/2- 

Translated into a condition on p, this ensures that H is self-adjoint. The 
corresponding solutions are 

^ 1 (-{-If eyiY>{\eE{^ - x^) + iknx)\ 

^ y exp(—^eE(=^ — x^) — i/cnx) j 

where kn = {n + ^jvr/L for n G Z. Before using them to construct a two- 
point function, let us pause to consider these solutions. Obviously, we have 
\Pnf‘{x) = SO the charge density of these modes is constant and does 
not depend on n. Hence, applying the mode sum formula and always 

^°Note the change of the sign at the two boundaries. It corresponds to contracting 7 
with the inward normal. 
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matching a positive and a negative mode, one obtains a vanishing vacuum 
polarization (filling positive modes or leaving holes in the negative modes 
would lead to a constant charge density). This seems to be the reason why in 
|19j . where the mode sum formula is employed, it is claimed that the vacuum 
polarization vanishes for massless fermions in 1 +1 spacetime dimension. As 
we will see, this is not correct. 

Setting = 'y^cpn, we may define the two-point function as in (31), (32), 
where the integrals are replaced by a sum over n, and the conditions on the 
integration region by kn ^ 0. To evaluate the current or the stress-energy 
tensor, we have to trace over the indices A, B with a 7 matrix. As these are 
off-diagonal, it suffices to compute the off-diagonal terms. We compute 

io{^pHx)My)) = 2 T/(f 

u{^Hx)My)) = 2x/(f -y°-x^ + y^- 

uj{My)i^\x)) = + y^-x^ + y^- 

uj{My)^\x)) = ^fili-x^ + y° + x^-y^- 

where 

f{^) = 

Hence, 


= — - — + 0{z^) 
2sinz/2 2 : 24 


uj{i> {x)il}2[y)) = 


2e2 


(x^ - y^)‘^{x^ -h 2 /^)^ 


27r 


ITT 


j;U _ yU _|_ 3,1 _ yi _ 


w(V'^(x)V'i(y)) = 


48L2' 


(34) 


2e2 


(x^ - y^)‘^{x^ -h 2 /^)^ 


27r 

iir 

" 48L2 


x^J — y» — x^ + y^ — ie 


(35) 


and analogously for the two-point functions (jj{'ipA{y)'ip^(x))- It is then clear 
that Schwinger’s prescription, with a coinciding point limit from the time 
direction, yields 

p(x) = ^e^Ex^, (36) 

which has the natural interpretation of negative charges accumulating where 
the potential is positive, and vice versa (note that the electrostatic potential 
is — Aq). However, omitting the parallel transport, one gets a vanishing vac¬ 
uum polarization p{x). This again exemplifies the importance of including 
the parallel transport in (29). 

The supplementary electromagnetic field produced by (36) via back- 
reaction is 

^br = ^e^E [(xi)2 - (L/2)2] . (37) 
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For the model to be physically viable, this should be small compared to the 
external field E, which leads to the requirement eL <C 1. 

In order to apply Dirac’s method, let us now compute the parametrix. 
In 1 + 1 spacetime dimension, the distributions T 2 k are, instead of (19), (20) 
given by 


T 2 =^(x,x') = 


-1 


22^-i7r(A: - 1)!^ A2 


.'^2 

±£ 


for fc > 1. Recalling (21) and (0 - ([g, we thus see that the parametrix 
is given by 


H^{x,y) = 


-i Vo{x,y)-f^{x - yY 


27r 


{x - y)±e 


up to terms that vanish as (x — yY log(a: — yY for a coinciding point limit 
from time-like or space-like direction. For the off-diagonal components, this 
means 


H^{x,yY2 = 


—il — ^eE(x^ + y^)(x^ — y^) — \{eEY{x^ + yY‘^{x'^ — y^Y 


2tt 


3;0 _ yO _|_ 3.1 _ yi ^ ie 


, 2 _ 1 “ \eE{x^ + - y°) - \{eEY{x^ + 


H^{x,yYi = 


27r 


— + ze 


again up to terms that vanish stronger than linearly in the coinciding point 
limit. Subtracting these from the two-point functions (34), (35), we obtain 
the remainder terms, cf. 


R{x, yY2 = i^eEix^ + y^ - Ye^^E\x^ + yY\x^ - y° - + yY 

-^{x^-y^ + x^-yY. 

E{x, yYi = ^eF;(x^ + yY - j^^e^E^ix^ + yY'^{x^ - y° + x^ - yY 
-^{x°-y^-x^ + yY, 

again up to terms of higher order. Taking the coinciding point limit and 
contracting with 7 ^^, one finds indeed a vanishing current Y the vacuum 
polarization (36). Note that filling supplementary positive energy modes 
(or dropping negative energy modes) just adds a constant charge density, 
which is independent of the energy level. The state that we have chosen is 
symmetric in the sense that the total charge density vanishes. In particular, 
this shows that the above results are not changed in a thermal state. 

Contrary to Schwinger’s prescription, Dirac’s method also allows for a 
computation of the stress-energy tensor. Its classical expression is given by 

5 ~ - VxYi^y) ■ 
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In particular, 


Too = I {'ip'yiVlip - Vi^7iV’) , 

^11 = I {'pjo^o'ip - Vo^ToV’) , 

Toi = I ('iA7iVoV’ + V’7oViV’ - Vo'tp'jiip - Vi'07oV’) • 


Evaluation in our state yields 


T^uix) = - 



/l 0\ 
^ y Vo V 


( 38 ) 


There are two contributions, one from the Casimir energy, and one from 
the interaction with the potential. The former coincides with the energy 
derived in [48| . The energy corresponding to the second term is — ^ times 
the electrostatic energy density of the charge density p in the potential Ao- 
Note that the integral over this energy density does not vanish]^ This shows 
that the popular mode sum expression, c.f., for example, the textbook [7j, 


^-tot 


1 

2 



is not correct, as this would yield Etot = 0. What may be more relevant than 
the energy density is the pressure at the boundaries, for which the electro¬ 
magnetic contribution is positive. This may be interpreted as a weakening 
of the attractive force of the two boundaries due to the shielding by vacuum 
polarization. Now there is the puzzling fact that we do not have 


dEtot 

dL 


P\±L/2- 


Instead, from (38) one obtains 


dEtot 

dL 


P\±L/2 — 


TT 

24L2 


2tt 


{L/2f, 


TT 

24L2 


2tt 


{L/2f. 


(39) 


The resolution of the puzzle is that one has to take into account the cor¬ 
rection to the electromagnetic stress-energy tensor due to the modification 


L —)• L -|- Lbr, cf. (37). At O(e^), this leads to an additional term 


Til = ^ [{x'f - (i/2)=] 

^^There is an ambiguity in the definition of the stress-energy tensor which amounts to 
adding multiples of the stress-energy tensor of the external fields, which is independent of 
. However, this ambiguity has to be fixed once and for all, so that it is not possible to 
achieve, for example, a vanishing total energy for all length L. 
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in the stress-energy tensor. Its contribution to the pressure at the boundary 
vanishes, but its contribution to the total energy is 


TTibr _ 

-^tot — 


2^ 2,^ ' ’ 


With the inclusion of this term, the relation (39) is indeed fulfilled. 


4 Static potentials in 3+1 dimensions 


In the previous section, we discussed an example where the two-point func¬ 
tion can be computed explicitly. In general, this will not be the case. Hence, 
some suitable approximation scheme is needed. This is greatly simplified if 
the external field possesses some symmetries. In the following, we consider 
the case of a static potential, i.e., = 0 and Aq is time-independent. 


For such a static potential, a ground state can be defined as in (31), 


(32). The practical problem is that in a numerical computation of the cor¬ 


responding vacuum polarization, the finite remainder of the coinciding point 
limit is difficult to extract. The idea will be to write the parametrix in a 


form similar to (31), so that one may integrate over the difference of the 


integrands. This gives an integral which is much better behaved. 

Due to the static potential, we are mainly interested in the vacuum 
expectation value of the charge density]^ Hence, let us consider the two- 
point function contracted with 7 *^ for a separation in time direction: 


Jo J E]^=E 


poo 

/ f{E,x)e-^^^dE. 

Jo 


(40) 


To arrive at the second line, we performed the integration and summation 
over momenta and spins for energy E (in the case of a radially symmetric 
potential, this will be a summation over angular momenta and spins). Now 
consider the parametrix contracted with 7 *^, as one performs the coinciding 
point limit from the time direction. Using the results of Section 2.1 
obtain 


we 




—2iVo(x + teo,x) i 2 ^o{x + te^,x) 

-— - 77 ? - 

vr^ E — ieE 2Tr^ t — ie 


e 

127r2 


d^Eox 



-)- iet 
A2 


+ 0{t). 


Up to the parallel transport, the first two terms also appear in the coinciding 
point limit of the contraction of the free two-point function, i.e., in the 

seems plausible that there are no rotational currents in the ground state, but we 
are not aware of a proof. 
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absence of external potentials. Denoting the sum of the first two terms by 
Ho{x,t), we thus have 


poo 

/ fQ{E,x)e-^^^dE 

Jo 


I 


ME - eAo{x),x)e-^^^dE 


eAo{x) 


where /o is the analog of / defined in (40) in the absence of external poten¬ 
tials. 

From the above considerations, it follows that we may write the expec¬ 
tation value of the charge density as 


p(x) = e lim 
t^+o 


[X[0,oo)(^)/(-E^,®) 

-X[eAo{x),oc){E)ME - eAo{x),x)] e~^^^dE 


(41) 




The second term in the integral has the natural interpretation of subtracting 
the vacuum contribution at the adjusted local energy level. Note that if 
the state w should be defined by a different choice of zero-point energy, 
one simply has to adjust the characteristic function multiplying f{E,x) 
accordingly. Also note that in some cases it may be advantageous to work 
with the two-point function for instead of the one for ip'll;. Then the 
sign in ED is changed, and the domains of integration are from —oo to 0 
(respectively eAo(^)). Finally, we note that in the absence of an external 
charge density, when the logarithmic divergence is absent, it may indeed be 
possible to perform the limit t —)• 0 in the integrand and have a numerically 
stable integral. However, this integral will in general still not be absolutely 
convergent, but oscillatory. 

We exemplify the above considerations with the Coulomb potential. To 
linear order in the external potential, the vacuum polarization was calcu¬ 
lated, using Dirac’s method, by Uehling M- Higher order corrections were 
later computed, using Schwinger’s prescription, by Wichmann and Kroll 
m, cf. also the review [3] and references therein. However, their calcula¬ 


tion did not take the parallel transport in (29) into account, which led to a 


divergent term at first order in the external potential, corresponding to the 


first term on the r.h.s. of (30). Being of the same order in the external field, 
this term was then identified with the Uehling term. Also at third order 
in the external potential some regularization is needed in the calculation of 
Wichmann and Kroll. Their procedure seems to boil down to also subtract 


the second, hnite, term on the r.h.s. of (30). So in the end, gauge invariance 


is restored, albeit in a rather intransparent way. 
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Figure 1: The vacuum polarization in a Coulomb potential (dots), computed 
using (41), as a function of the radius (in units of the electron Compton wave¬ 
length) for aZ = 0.1. As a comparison, the vacuum polarization according 
to Uehling is also shown. 


In order to illustrate the usefulness of the expression (41), we use it to 
compute the vacuum polarization in the Coulomb potential. The explicit 
form of the modes of the continuous spectrum can be found for example 
in [39]. As the external charge density vanishes except at the origin, the 
second term on the r.h.s. of (41) is absent, and the limit t —)• 0 may be 
performed inside the integral. We consider the state obtained by Filing up 
the negative continuum, i.e., by considering a completely ionized nucleus. 
Figure shows the charge density p{r) for aZ = 0.1, as a function of the 
radius. We see a nice agreement with the result of Uehling m for the 
vacuum polarization linearized in the external field. These numerical results 
mainly serve as an illustration that the above method can be practically 
applied. The concrete numerical implementation can certainly be improved 
considerably. In particular it should be possible to compute the higher 
order corrections to Uehling’s result. In a separate publication [30|, a time- 
dependent homogeneous electric held is treated in an analogous way, with 
agreement to well-established results both in the perturbative and the non- 
perturbative regime. 
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